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ABSTRACT 

Narrowband  filters  are  inherently 
of  high  order,  so  they  require  a high 
computation  rate  if  implemented  directly. 
Multirate  filters  are  capable  of  effec- 
tively approximating  narrowband  charac- 
teristics. Such  an  implementation  may 
require  a much  reduced  rate  of  computa- 
tion. This  paper  derives  an  expression 
for  the  system  error  of  this  approach. 
This  error  consists  of  two  components, 
one  due  to  the  filter  approximation  and 
the  other  due  to  aliasing.  The  role  of 
aliasing  error  in  the  design  of  multi- 
rate filters  is  then  considered.  Ex- 
amples are  presented. 


I . INTRODUCTION 

Multirate  filters,  those  filters 
composed  of  cascaded  decimators  and  in- 
terpolators, have  been  shown  to  be  useful 
in  implementing  narrowband  lowpass  fil- 
ters. Bellanger  et  al.  in  [1]  used  multi- 
rate filters  composed  of  half  band  deci- 
mators and  interpolators  for  lowpass  fil- 
tering and  they  demonstrated  that  compu- 
tational savings  could  be  realized.  Rab- 
iner  and  Crochiere  in  [2]  and  [3]  found 
multirate  lowpass  implementations  which 
minimized  either  multiplication  rate  or 
coefficient  storage.  Both  of  these  works 
assumed  that  the  aliasing  error  was  ne'g- 
ligible. 

In  this  paper,  we  calculate  the  sys- 
tem error  of  multirate  implementation  of 
lowpass  filters.  This  error  consists  of 
two  compoft^ts,  one  due  to  filter  approxi- 
mation and  tl^e  other  due  to  aliasing.  By 
using  this  re^lt,  design  rules  can  be 
formulated  whic>i  limit  the  aliasing  error 
in  these  filters.  Examples  of  multirate 
filters  so  designed  are  given. 

^This  research  is  sponsored  by  the  Air 
Force  Office  of  Scientific  Research,  Air 
Force  Systems  Command,  USAF,  under  Grant 
No.  AFOSR-76-3083.  \ 


II.  THE  SYSTEM  ERROR  FOR  AN 
N-STAGE  MULTIRATE  FILTEk 

Let  £x^)  be  the  input  to  and  (y**} 

the  output  of  a filter  with  tremsfer  func- 
tion F^(z),  as  illustrated  by  the  upper 

signal  path  of  Fig.  1.  An  N-Stage  multi- 
rate implementation  of  this  filtering, 
shown  as  the  lower  signal  path  of  Fig.  1, 
is  composed  of  N generalized  decimators 
followed  by  N generalized  interpolators. 

Its  output  is  ly  ] which  is  to  be  close 

to  ty^). 

A generalized  decimator  consists  of 
a filter  tt^H^(z)  followed  by  a decrease 

in  the  sampling  rate  by  a factor  of  D.. 

It  reduces  to  a simple  decimator  if  Hj^  is 

a lowpass  filter.  A generalized  interpo- 
lator consists  of  an  increase  in  the  sam- 
pling rate  by  a factor  of  D^  followed  by 

the  filter  g.G. (z) . 

11 

Let  H^  and  be  linear  phase  FIR 
filters  of  orders  Njj^  and  N^^  respectively. 
The  H^  are  symmetric,  and  only  one  output 
of  H^  is  required  for  every  D.  inputs. 

Only  one  of  each  D^  inputs  to  filter  G^  is 

non  zero.  By  taking  advantage  of  these 
properties,'  the  raultirate  filter  requires 
N i 

only  E (N_.  + N ./2)  / 1 f D.  multiplica- 

i=l  j=l  J 

tions  per  output,  if  this  is  less  than  the 
multiplication  rate  of  a comparable  direct 
implementation,  a computational  saving  is 
realized  by  using  the  multirate  structure. 

In  the  lower  path  of  Fig.  1,  the  sam- 
pling rate  is  increased  and  decreased  ac- 
cording to  the  equations : 
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(1) 


2i  _ 2i-l 

^n.  ' ’'o.n.  +9. 
1 11  1 


4N-2i+l 


n.  1 =D.  n . +9  . 
1-1  1 1 • 1 


, otherwise 


for  i=l,2 N where  6^^  is  an  integer 

valued  random  variable  taking  on  the  values 
0, 1 , . . . , Dj^-1  with  equal  probability. 

Q For  a wide  sense  stationary  input 
{x^) , with  power  spectral  density  ♦ [exp  (jiu)  ] , 
it"  is  shown  in  [4]  that  the  mean  square 
error  is  given  by 

E{  (y„-y^)^) 

TT  N 2 

= Tl  TT  H.  [g.  (ui)]G^  [g^^  («)]  - Fjtgj(iu)]| 

-Tf  i=l 

X ♦ [g,  (uj ) ] diu/2TT 
n N , 

+ Jlr-  TTg.  [g.  (ui)]H  [f.  (U)))r#  [f^(m)]<3m/2n 
i=l  (3) 

for  Tfa.e.  = TT'5-  where 

i=l  ^ ^ i=l  ^ 

g^(n))  = exp(jmDj^. . .D^_3)  . (4) 

i-1  N k 

f.  (hi)  = cxp(jiu  TTD.+j2f»  C n./  T~fD.  ) , 

j=l  ^ k=i  ^ j=i  3 (5) 

and  the  product  [Dj^ . . .D^_  is  defined  to 

be  1 for  i-1.  The  symbol  represents 

N summations?  one  for  each  n^.  The  range 

of  summation  of  n^  is  from  0 to  D^-1,  but 

the  term  where  all  n.  equal  zero  is  ex- 
cluded from  the  summation. 

The  first  term  in  (3)  is  the  error 
due  to  approximation  of  F,  [expCjw)]  by 
N 

TT  [gj^  (hi)  (q^  (ui)  1 . We  call  this  the 

filter  approximation  error  and  denote  it 
by  e . The  second  term  in  (3)  is  the  er- 
ror ^ due  to  aliasing  generated  by  the 
sampling  rate  changes.  We  call  this  the 
aliasing  error  and  denote  it  by 

For  acceptable  filter  approximation, 
we  require  that  the  magnitude  of 
N 

TT H.  [g.  (h))]G.  [g.  (<u)]  be  within  6 of 
1 ^ ^ ^ P 

Fj[exp(jui)]  in  its  passband  and  within 

6 of  F,  [exp(jm)]  in  its  stopband.  For 
& a 

acceptable  aliasing  error  we  require  that 
e^  * E{  (x°)^)/S  (6) 

for  some  specified  S.  Since  the  aliasing 


is  correlated  to  the  signal  it  is  more  ob- 
jectionable than  noise  which  is  not.  Thus 
a reasonable  S would  make  e,  comparable  to, 
but  somewhat  less  than  the  roundoff  noise. 

By  picking 

“•sHi  = “sGi  = 2"/°i  - °l-”'^i-l“’s  ' 

“sHN  = '^GN  = °l'-'Vl‘"s  ' 
an  approximation  Cor  e^^  can  be  made  which 

is  valid  for  a reasonably  smooth  ♦[exp(jn))]. 
This  approximation  allows  us  to  control  the 
aliasing  error  by  simple  choices  of  the 
stopband  ripples  of  H^[exp(j'u)]  and 

G^[exp(jw)].  'The  procedure  of  Crochiere 

and  Rabiner  in  [5]  is  us^d  to  pick  the  D. 
subject  to  the  restrictions  of  (7)  and  ^ 

(8) . The  rest  of  the  design  parameters  of 
the  H.  (exp(jii))]  and  the  G.  [exp(j(u)]  are 

then  picked  to  satisfy  the  filter  approxi- 
mation requirement.  The  details  can  be 
found  in  [4] . 

III.  EXAMPLES 

We  present  here  three  examples  of 
multirate  filter  design. 

Example  1 . The  specifications  are: 

W =.05tt,  6 =.01,  u)_=.ltT,  6 =.001,  S=10^. 
p p s s 

By  using  the  design  program  of  McClellan 
et  al.  [6],  a filter  of  order  of  110  was 
found  to  satisfy  the  specifications.  The 
design  parameters  for  a one  stage  multi- 
rate filter  are: 

i»pH=.05n,  6pH=.0n5,  mgjj=.ln,  6gjj=.0046  ; 

III  _=.05n,6  =, 005,111  _=.!«, 6 _=.0015;D=10. 

pG  pG  SG  sG 

Filters  were  designed  to  these  specifica- 
tions. The  orders  of  H and  G are  96  and 
•110  respectively.  The  optimization  pro- 
cedure of  Crochiere  and  Rabiner  in  [5]  se- 
lects as  optimal  Dj^=5,D2=2.  For  this  two 

stage  implementation,  our  filter  design 
parameters  are 

“pHl^ • ' ®pHl“ • ^ ' '"sHl=  • ' '’sHl“ • ' 

'“pGl°-”^"'W°-°"^^'“’sGl=-30"«6sGl=-°"^®' 

‘''pH2”  -25".  • 002 5 . »g jj2  = • son , 6g = • 0046 ; 

•pG2=-25n,  6pH2=-0025,ii;sG2  = - son,  6sg2“*0O18. 

Filters  were  designed  to  meet  these  para- 
meters using  [6].  The  orders  of  H, ,G, ,H, 
and  G2  were  20,22,23,25.  112 

Figure  2 illustrates 
N 

7T  H.  [g.  (w)]G.  [g.  (ai.)]  for  the  one  stage 
i=l  1 1 

implementation.  Figure  3 illustrates  it 
for  the  optimal  two  stage  implementation. 
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Table  1 gives  a comparison  of  the  multi- 
plication rates  of  the  various  implementa- 
tions . 

Example  2 . The  specifications  are 

tu  =.0833tt,6  =.31.ii)  = 1667«,6  =.155,S=10^. 
p p s s 

The  Crochiere  and  Rabiner  procedure  gen- 
erates Dj^=3,  optimal  choices 

of  the  Table  1 compares  computation 

rates  for  the  direct,  one  stage,  and  opti- 
mal implementations. 

Example  3.  The  specifications  are 

(U  =.0833n,6  =.31,tt)  =.1667n,6  =.155,S=10^. 
p p s s 

Table  1 compares  the  multiplication  rates 
for  the  direct,, one  stage  and  two  stage 
implementations . 

IV,  CONCLUSIONS 

Aliasing  error  arises  in  a multirate 
filter  as  a result  of  sampling  rate 
changes . We  have  calculated  the  aliasing 
error  in  a multirate  filter,  and  designed 
multirate  filters  which  impose  a bound  on 
aliasing  error  while  meeting  the  filter 
approximation  specifications.  Examples 
2 and  3 demonstrate  that  this  multirate 
filter  may  not  have  a computational  ad- 
vantage over  the  direct  implementation. 

In  many  cases,  however,  the  multirate 
filter  will  have  a computational  advantage 
as  it  does  in  Example  1. 
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Fig.  1.  N- STAGE  MULTIRATE  IMPLEMENTATION 

(top  path:  desired  filtering 
bottom  path:  multirate  implementation) 


Fig.  2.  FREQUENCY  RESPONSE  FOR  THE  Fig.  3.  FREQUENCY  RESPONSE  FOR  THE 

ONE- STAGE  IMPLEMENTATION  OF  EXAMPLE  1 TWO- STAGE  IMPLEMENTATION  OF  EXAMPLE  1 
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20.  AbrjtruOi 


>fone  due  to  the  filter  approximation  and  the  other  due  to 
/ aliasing.  The  role  of  aliasing  error  in  the  design  of  multi- 
rate filters  is  then  considered.  Examples  are  presented. 
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